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Abstract 



We revisit and extend previous work on neutrino mediated long range forces in a 
background at finite temperature. For Dirac neutrinos, we correct existing results. 
We also give new results concerning spin-independent as well as spin-dependent 
long range forces associated to Majorana neutrinos. An interesting outcome of the 
investigation is that, for both types of neutrinos whether massless or not, the effect 
of the relic neutrino heat bath is to convert those forces into attractive ones in the 
supra-millimeter scale while they stay repulsive within the sub-millimeter scale. 



Neutrinos mediate long-range forces between macroscopic bodies |]T|, 0, 0, 0], p. In- 
deed double neutrino exchange among matter fermions generates spin-independent forces that 
extend coherently over macroscopic distances. The effect, however, is extremely weak, much 
too tiny to be experimentally detected with present day technology. Compared to their grav- 
itational pull, the force between two nucleons 1 cm apart is about 10^^^ times weaker. Not 
only their coupling strength is very small but also their decay with distance is fast. Indeed 
the potential drops as so that the effects die off correspondingly. Phenomenological sur- 
veys on forces with this particular distance behaviour have been conducted in the literature 
(see e.g.0) over the whole span of distances from astronomical down to the micron scale. If 
at all, this forces will induce physical effects in the sub-millimeter (but macroscopic) end of 
the distance scale. Perhaps an exception to this is the case of a system with high density of 
matter such as the core of a neutron star where collective effects may show up |^, [^. 

In a neutrino populated medium, such as the cosmic neutrino background or the hot core 
of a supernova, the helicity flip produced by single neutrino exchange can be balanced by the 
neutrinos in the medium and, as a consequence, a spin-independent interaction takes place 
that leads to a coherent effect over many particles in bulk matter. 

The neutrino long-range forces in the presence of a neutrino thermal bath have been 
explored in reference p] in the Dirac neutrino case. The long range forces mediated by 
Majorana neutrinos, on the other hand, have been studied only in the zero temperature case 
[^. Here we wish to extend the nonzero temperature results to the Majorana case. Because 
the distinction between Dirac and Majorana neutrinos is superfluous for massless neutrinos, 
we shall consider the general m 7^ case. 

We shall adopt the notation in |^ and write. 



where T(Q) is the nucleon-nucleon elastic scattering amplitude (figure 1) in the static limit, 
i.e. momentum transfer Q ^ (0, Q), where matter is supposed to be at rest in the microwave 
background radiation (MWBR) frame. It can be cast in the form 




(1) 



T(g) = -2iGU9v,-2gASng'y,-2g'^Srif 
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with 
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Figure 1: Fig. 1. Lowest order Feynman diagram for two neutrino exchange in the four 
fermion effective theory. 

The operator O is the left-handed projector |(1 — 75) for Dirac neutrinos and -^75 for 
Majorana neutrinos. The temperature dependent propagator St has the exphcit form 

ST{k) = {f + m) - + i^y^ + 2m5{k'^ - m2)(^(A;°)n+ + ^(-A;°)n_)] (4) 

where n_|_ and n_ are Fermi-Dirac distribution functions for particle and antiparticle, re- 
spectively. As discussed in [Q, figure 1 evaluated with this propagator taken together with 
the usual Feynman rules is sufficient to calculate the potential. In equation (2), gy^A are 
composition-dependent weak vector and axial-vector couplings. We focus first on the spin- 
independent potential, that is the gvQy component of equation (2). 

Use of the first piece in equation (4) gives the zero temperature vacuum results 0, P, 

= ^^£^^3(2mr) (5) 

and 

VMajoranair) = — ^^-^-^-^ ( 2mr ) (6) 



in terms of the modified Bessel functions K2. 
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At very large distances (mr ^ 1), i.e. much larger than the Compton wavelength of the 
neutrino, these potentials exhibit the asymptotic behaviour 

Vo,™.(r) (^)''%— (7) 

and ^ ^ 

. (r) ~ ^f9v9v ( ^\ -2mr /oN 

Of course, both potentials (equations (5) and (6)) coincide when m = 0. They give the well 
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known Feinberg and Sucher result: 



In a neutrino background, a contribution to the long range force can arise because a neutrino 
in the thermal bath may be excited and de-excited back to its original state in the course of 
the double scattering process. This effect is described by the crossed terms contained in 7^,^ 
that involve the thermal piece of one neutrino propagator along with the vacuum piece of the 
other neutrino propagator. This thermal component of the tensor I^i, can be written as 



X 



. (27r)4 



{k + Qy — m? + ie {k — QY — m? + ie 



(10) 



in the Dirac case, and 



n 



X 



d^k 

Tr [-"{p + (^ + ;n)Y{P - ni 



{k — QY — rn^ + ie 



(A; + Qf -m^ + ie 

for Majorana neutrinos, where in this latter case wc put n+ = n_ = n since the chemical 
potential vanishes. Note that in (10) there is no component proportional to e^'^°'^ since after 
the integration over k the only four- vector available is Q". As a result there will be no parity 
violating potentials. 

Far from degeneracy (i.e. for chemical potential // <S T), as is probably the case for 
cosmological neutrinos, we can consider the neutrinos to be Boltzmann distributed, that is 
we take 

n± = cxp[(±/x-|A;°|)/T]. (12) 

With this approximation, the integrations involved in the calculation of potentials can be 
easily done by conveniently choosing the order in which they are performed. The results can 
be expressed again in terms of Bessel functions and are as follows: 



yDirac, 



r) 



Gj,m^gvgv 



TT^r 



cosh {iJ>/T) 



P 



P 



and 



yMajc 



(13) 



(14) 
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where we have defined 



m 



p=-Vl + (2rT)2. (15) 



For massless neutrinos (and = 0) both potentials collapse to 

r = § 1 16 

which is the result given in reference [^. Because the neutrino background temperature is 
T ~ (1.2mm)^^, we see that for distances much larger than 1mm (i.e. rT ^ 1) the potential 
in equation (16) reads 

Gf9v9 v 

When added to the vacuum result (9), the total potential is 



Vrir) ^ (17) 



V.Ar) - (18) 

that is, in the presence of the cosmic neutrino background the original Feinberg-Sucher force 
switches sign, i.e. a repulsive force turns into an attractive one. On the other hand, well within 
the sub-millimeter domain (rT -C 1), the temperature dependent potential (16) behaves as 
follows 

V^{r) ^ (19) 

which is negligible compared to the vacuum contribution in equation (9). 

In the general m ^ case, we shall study the Dirac and Majorana potentials, equations 
(13) and (14) respectively, in various physically interesting limits. Consider first the cases 
where r ^ m^^ ^ T^^ or m^^ ^ r ^ T^^. Performing the relevant expansions of the Bessel 
functions in (13) and (14) leads to 

5/2 / 

(r) ^ T^" cosh (/i/T)e-™/^ (20) 

and 

Hence, thermal effects are exponentially damped in both distance domains. 

A different behaviour is obtained for distances much larger than any inverse energy scale 
in the problem, i.e. for r ^ T^^ ^ m^^ or r ^ m~^ ^ T"^. Indeed, now we have 

^?--(r) ^ -9kx9v f I^y/' ^^^^ (p/T)e-2"^^ (22) 
4 Vvrr / 

and 

y Majorana ^ _Gp9v9y_ f \ ^_2mr j^23) 



Both expressions exhibit the characteristic Yukawa exponential damping associated to two- 
particle exchange. These results when added to their vacuum counterparts, equations (7) and 
(8), produce the inversion phenomenon already noticed in the massless case. At asymptotically 
large distances the resulting potential is equal in strength as it would be in vacuum but, 
contrary to what happens in vacuum, it is attractive instead. 

There is no exponential suppression only when r <^ T^^ ^ or ^ r ^ m~^, 
where one essentially recovers the massless cases, equations (19) or (17), respectively. Indeed, 
for Majorana neutrinos one gets these equations as they stand, and for Dirac neutrinos both 
equations should be multiplied by the factor cosh (yu/T) for non-zero chemical potential. 

Let us note that the results given in equations (20) and (22) for the m 7^ Dirac case 
disagree with the corresponding results given in reference [^. Indeed, their formulae do not 
show the Boltzmann or Yukawa suppression factors that enter the asymptotic expansions of 
the Bessel functions and which are bound to be there on physical grounds. For the sake of 
an explicit comparison we provide the reader, in the appendix at the end of the paper, with 
some details of the calculation. 

Up to this point all calculations refer to spin-independent potentials, those that can coher- 
ently add over macroscopic samples of unpolarized matter. Let us, for the sake of complete- 
ness, consider briefly the question of potentials that depend on spin. Now we should focus 
on the spatial indices of the tensor J^,^ appearing in the scattering amplitude in equation (2). 
The Fourier transformation (1) is in this case somewhat more involved than before because 
the amplitude will depend on the components of the 3-momentum transfer. Nevertheless they 
can be easily performed and we get: 



spin/ 



AGlm^gAg'A 



[S-S') F(r) + 2 ^^"'^f "^^ G(r) 



cosh(/i/T) (24) 



where 



and 



F(r).a^ + 2^-8.nV^ 



G{r) 



P P' 
^ ^K2{p) 



. 2 2^Mp) 

4m r — 



(25) 



(26) 



with a = 1 for Dirac neutrinos and a = 2 for Majorana neutrinos and = 3/4. Of course, 
in the Majorana case we must put fi = 0. 



Both cases above, i.e. Dirac and Majorana, lead to the potential 



vr3(i + 4r2T2)3r 



S ■ S')(l - Ur'T') + (^•^)(S'-^) (7 + i2r2T2) 



(27) 



5 



for m = and ^ = 0. This result, eq. (27), should be then added to the vacuum result [|I] 

Glgyg'y r (ST)(S--r) _ _ 



{2i 



The various regimes explored before can be studied also for the spin- dependent forces. The 
discussion involves the various asymptotic forms of the same modified Bessel functions and 
will lead to the same exponential damping whenever the temperature or the mass is the 
relevant energy parameter. Since these forces will be even more difficult to detect than the 
spin-independent ones, for they do not add up coherently in bulk matter, we do not bother 
here to display the explicit form for the different limits. 

We end this paper with a short summary. Double neutrino exchange mediates (extremely 
feeble) long range forces. In vacuum these forces have been known (at least for Dirac neutri- 
nos) for quite some time. Recently, it has been realised that a neutrino background will also 
induce long range interactions among bulk matter. The results were given for Dirac neutrinos. 
We have extended the work of Horowitz and Pantaleone to include the case of Majorana 
neutrinos and, furthermore, we have derived the exact form of the potentials in either case, 
i.e. Dirac and Majorana, and explored physically relevant distance and energy scales. In so 
doing we have found important discrepancies with previous work (m 7^ 0, Dirac case). Since 
matter is embedded in the cosmic neutrino background, a consequence of our analysis is that 
the forces are repulsive in the sub-millimeter scale and attractive for distances well beyond 
1mm for any kind of neutrino (massless or not). In fact, on the small scale the vacuum result 
(Feinberg and Sucher) dominates whereas on the larger scale the relic neutrino background 
is responsible for the dominant effect. This means that by experimentally detecting (ad- 
mittedly a highly improbable event for laboratory experiments) such forces in both different 
regimes one would, not only establish these neutrino interactions, but one would in addition 
detect the relic neutrino background. Actually it is the neutrino background temperature 
(T^^ ~ 1.2mm) which sets this 1mm distance scale. Incidentally, the sub-millimeter scale 
has been subject recently of renewed theoretical as well as experimental interest [0. For an 
experimental point of view of the actual possible detection of very weak long range forces we 
refer the reader to references ||1CI|| . 
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Appendix 



We give here some details of the calculation for the massive Dirac case. All integrals used can 



be found in |rT|. Starting with eq. (1) we have (in the formulae below Q = |Q|): 



Dirac ( 
T \ 



dQQIooiQ) sin(Qr) 



cosh {jJi/T) 



dkk^ 



X 



\^W+m? 



exp (^—Vk'^ + m? /T^ 
dQ 



Q2 _ {2kzY 

Performing first the integration over Q, then over we obtain: 



Dirac i 



r) 



cosh (/i/T) 



— (1 + (mr)^) /(m, r, T) + r 



dl{m, r, T) 
dr 



where 



/(m, r, T) 



dkk 



exp + m? /T^ sin(2A;r 



\/k'^ + m? 
l + (2rt)2 



m 



Inserting (A. 3) in (A. 2) gives eq.(13). 

For the spin-dependent part we decompose 



T(Q) = (S • Q)(S' ■ Q) h{Q) + (S ■ S') UQ) 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



where the functions ti{Q) and t2{Q) depend only on |Q| and we perform first the angular 
Q integration in eq.(l). The rest of the calculation goes along similar lines as in the spin- 
independent part above. 
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